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SUMMARY 

A  computer  study  is  described  of  the  effect  on  I.L.S.  glide  path 
performance  of  the  addition  of  various  damping  terns  to  the  glide  path  control 
law  of  a  typical  autopilot.  The  control  laws  used  are  analysed  by  a  rwnfcor 
of  techniques  arri  the  results  achieved  by  the  different  techniques  are  compared. 
It  is  shown  that  the  performance  can  bo  improved  by  the  addition  of  either  a 
vortical  velocity  term  coupled  with  an  acceleration  term,  or  a  vortical  volocily 
term  coupled  with  a  pitch  rate  term  in  the  control  law. 
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1  INTRODUCTION 

One  of  the  more  important  performance  parameters  of  an  autonatio  landH* 
system  is  the  longitudinal  position  of  the  touchdown  point  of  the  aircraft 
along  the  nrmay.  The  mean  position  and  the  scatter  of  the  touchdown  point 
should  be  such  that  the  chance  is  extremely  remote  of  the  aircraf  • 
larding  short  of  the  runway  threshold,  i.e.  undershooting,  or  floa  ng  ** 
the  runway  so  that  it  over-runs  the  far  end.  VOrilst  the  quality  of  the  flare- 
out  manoeuvre  can  affect  the  position  of  the  touchdown  point,  the  accuracy 
performance  achieve  on  the  I.L.S.  glide  path  Wiately  prior  to 
has  a  far  grater  effect.  The  positional  error  of  the  aircraft  fnom  the  I.L.S. 
glide  path  aid  the  rate  of  change  of  this  error,  which  is  a  measure  of  the 
error  in  the  aircraft's  rate  of  descent  fran  that  required  to  fly  down  the 
glide  path,  each  have  a  considerable  effect  on  the  point  of  touchdown  . 


.\n  assessment  of  the  recorded  results  obtained  from  full-scale  flight 
trials  has  shown  that  the  longitudinal  soatter  of  the  touchdown  point  of 
current  automatic  larding  systems  can  only  be  constrained  within  the  limit 
imposed  by  Civil  safety  standards  by  restricting  the  wind  conditions  in  which  the 
systems  are  used.  For  operational  use,  wider  wind  limits  are  desirsble  and  so 
an  analogue  computer  study  has  been  made  to  investigate  whether  improvements  can 
be  made  to  the  accuracy  of  performance  on  the  I.L.S.  glide  path.  This  Report 
describes  the  results  of  the  investigation. 


One  way  of  improving  beam  holding  in  the  presence  of  wind  disturbances 
is  to  increase  the  gain  of  the  coupling  between  the  aircraft  system  aid  the 
radio  beam.  However,  increase  of  gain  alone  causes  a  loss  of  overall  system 
stability  with  the  result  that  any  improvement  in  performance  by  this  manner 
is  extremely  limited,  in  faot  the  performance  is  more  likely  to  deteriorate. 
Hovvovor,  the  use  of  phase  advance  toms  can  increase  systom  stability  aid  ao 
allow  an  increase  of  coupling  gain,  and  the  addition  of  such  terms  to  the 
basic  glide  path  control  law  of  a  typical  autopilot  has  been  investigated  on 
a  ccmouter.  The  terns  used  have  ranged  from  the  time  derivative  of  the 
glide*  path  signal  to  combi  rations  of  terns  obtained  from  a  normal  accelerometer 

aid  pitch  attitude  gyros. 

In  this  Report,  emphasis  is  laid  on  frequency  response  and  statistical 
response  techniques  of  analysis,  since  difficulty  has  bean  encountered  in 
the  past  in  the  prediction  of  performance  from  time  histories  of  simulated 
approaches  in  the  presence  of  steady  or  slowly  varying  wind  conditions. 

The  criterion  against  which  performance  is  measured  is  the  performance 
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aohievod  using  the  simulation  of  the  u  modified  control  law  of  the  Mark  10B 
autopilot. 

Seotion  2  discusses  the  analysis  procedures  adopted,  together  with  the 
techniques  used  for  choice  of  gearings  and  the  difficulties  associated  with 
those  and  past  techniques.  Section  3  gives  a  statement  of  the  modifications 
investigated,  followed  by  a  detailed  diecussion  of  results  in  Section  4. 

The  modified  glide  path  control  equations  used  in  this  study  are  listed 
in  the  Appendix  with  the  longitudinal  equations  of  motion  of  the  aeroplane 
(Varsity)  and  th^  kinematic  and  wind  equations. 

2  ANALYSIS  PROCEDURE 

In  tho  past  there  have  been  marked  discrepancies  between  computer  studies 
of  control  on  the  I.L.S.  glide  path  and  the  results  of  actual  flight  trials. 

For  exaitple,  work  done  using  the  time  history  record  and  single  discrete  wind 
gusts  as  disturbances,  erroneously  predicted  a  possible  gain  factor  of  two  or 
three  times  the  highest  practical  flight  value.  For  some  time  these  discrep¬ 
ancies  were  blamed  on  inaccurate  or  ircomplote  simulation  of  the  problem,  but 
recent  work  has  indicated  that  inadequate  methods  of  assessing  the  performance 
achievod  on  the  computer  were  the  main  cause  of  the  differences.  One  difficulty 
lies  in  the  range  dependence  of  the  glide-path  signal.  The  contours  of 
constant  error  signal  are  straight  linos  radiating  from  the  transmitter,  so 
that  the  error  signal  received  by  the  aircraft  for  a  given  linear  displacement 
from  tho  glide-path  beam  is  greater  at  short  range  than  at  long  range  from  the 
transmitter.  Thus  the  system  gain  is  low  at  long  range,  giving  sluggish  control, 
but  high  at  short  range,  giving  a  tight  control.  In  this  situation  it  is 
inadequate  to  examine  the  time  history  of  an  approach  in  which  the  only  distur¬ 
bance  is  the  error  at  the  start,  since  by  the  time  the  system  reaches  instabi¬ 
lity  because  of  gain  increase,  there  will  be  no  disturbance  to  trigger  the 
instability.  Similarly  the  injection  of  a  step  disturbance  at  a  given  range 
is  inadequate  since  the  effect  of  this  disturbance  is  critical  upon  the  range 
at  which  it  is  injected,  and  it  is  not  clear  how  to  choose  a  range  that  will 
be  characteristic  of  the  whole  approach. 

It  is  clear  that  the  only  satisfactory  test  is  to  subject  the  system  to 
rardom  disturbances  throughout  a  series  of  approaches  (as  would  occur  in 
practice)  ani  to  examine  statistically  the  results  of  such  a  series.  Subsid¬ 
iary  tests  could  be  used  to  indicate  whether  the  beam  capture  manoeuvre  was 
satisfactory. 
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The  method  of  analysis  adopted  in  this  Report  falls  into  four  diatinot 
stages,  of  which  the  first  three  are  additional  to  previous  mottos 

(i)  choice  of  control  gearing  values  to  be  used  in  the  modified  control  laws 
by  a.  study  of  stability  margins; 

(ii)  frequence  responses  of  displacement  fran  the  I.L.S.  glide  path  cau"®* 
hy  three  variables  -  (a)  vertical  wind,  (b)  horizontal  wind  (c)  radio  no  , 

(iii)  statistical  results  from  sinulated  approaches  in  the  presence  of  rw.dom 
horizontal  wind  variations; 

(iv)  «  study  of  opprooch  mooumroo  in  ths  pro.unou  of  vsrious  dlsturksoooo, 
„oh  so  th.  loo.  join  imnoouvrs ;  «ho  tffoot  of  rot.  of  «~°0»t  «•*»  «<“• 
and  the  effect  of  wind  shears  both  from  head  and  tail. 

Pig. 10  shows  a  block  schematic  of  the  computer  arrangement. 

2.1  Choice  of  gearing  values 

Since  it  is  pointless  to  have  a  system  which  gives  good  performance  at 
long  range  but  is  unstable  at  close  range,  the  gearing  values  chosen  were 
those  which  gave  as  high  an  overall  gearing  as  was  feasible  at  close  range, 
since  high  gearings  produce  tighter  control  in  the  presence  of  wind  turbulence. 
Tho  general  method  for  choice  of  gearings  is  illustrated  in  the  following 
example.  With  a  general  control  equation  of  tho  form:- 


0.  =  -K. 


(» ♦  I  *  K10,  “) 


where  6  =  commanded  pitch  angle, 

c 

'3  =  glide  path  error  signal, 

H  =  height  of  aircraft, 

D  -  ft 

K_,  K,,  K...  are  control  gearings, 
y  6'  101 


assume  that  tho  required  values  are  those  of  K$  and  K^,  the  value  of 
being  invariant.  Then  K  is  fixed  at  a  particular  value  and  K5  varied 
until  the  recovery  of  the  system  from  a  step  disturbance  is  a  neutrally  domped 
oscillation.  The  procedure  is  repeated  at  other  values  of  and  the  value 

of  K,  for  neutral  stability  plotted  as  a  function  of  K^,  thus  giving  the 

neutral  stability  bourdary  in  tho  pl«».  Sirica  :nininum  BtaMlity 

occurs  at  tho  start  of  tho  ATTITUDE  phase  (100  ft)  when  tho  coupling  to  the 
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glide  path  in  terms  of  commanded  pitoh  attitude  per  foot  of  displacement  from 
the  glide-path  beam  is  a  maximum,  the  above  boundary  is  computed  with  a  fixed 
height  of  100  ft  programmed  into  the  computer. 

Stability  boundaries  similar  to  the  above  have  been  plotted  for  each  of 
the  different  control  equations  listed  in  the  Appendix,  and  the  boundaries  are  shown 
in  Figs.1  to  7.  It  was  found  that  in  the  case  of  DH  damping  alone  (Fig,  l)  the  curve 
reached  a  maximum,  but  for  most  of  the  other  modifications,  (Figs, 2  to  5)  the  curve 
was  a  monotonic  increasing  function  of  the  gearing  K1Q1  (or  K 1Q2,  etc)  up  to  a 
point  where  the  rate  of  pitch  command  exceeded  the  limit  of  3°/sec,  when  the 
neutrally  damped  oscillation  suddenly  degenerated  into  a  rapidly  divergent 
oscillation  at  the  slightest  provocation,  Tho  valuo  of  Kc  at  which  this  sudden 
chango  took  place,  was  not  consistent  for  a  given  value  of  K ^  (or  K^,  etc)  and 
this  area  of  uncertainty  is  shown  by  the  hatched  aroas  of  Figs.  2  to  5.  This  effect 
initially  caused  concern  as  to  tho  restrictions  it  would  apply  to  use  of  this 
method  of  choosing  gearing  values  and  an  inquiiy  was  made  as  to  its  nature. 

With  no  rate  limit  present  the  boundary  continued  anoothly  with  no 
discontinuities,  Addition  of  the  rate  limit  alone  caused  the  system  to  become 
divergent,  but  limiting  the  pitch  command  amplitude  to  its  standard  value  of 
3?°  contained  the  oscillation,  although  it  did  not  damp  it  out. 

While  the  possibility  of  such  limit  cycles  might  appear  unsatisfactory,  it 
was  found  that  with  tho  rate  and  amplitude  limits  both  present  on  tho  pitch 
coermond  and  tho  finally  cho"  m  goarings,  an  offset  of  ft  from  tho  glide-path 
at  a  height  of  100  ft  was  required  to  initiate  the  undamped  oscillation,  and 
at  greater  heights  the  required  initial  offset  was  greater  than  the  beam  width. 

This  was  considered  to  be  a  condition  most  unlikely  to  be  met  in  practice  since 
both  rate  and  amplitude  limits  arc  inserted  for  reasons  of  passenger  comfort 
and  safety. 

Flight  cxj)crionce  has  shown  that  the  optimum  value  of  the  terminal  gearing  i3 
about  half  that  value  which  gives  instability.  The  choice  was  therefore 
nominally  half  the  maximum  value  of  attainable  on  tte  stability  boundary 
plots.  Unfortunately,  due  to  the  effects  of  the  limit  cycles  referred  to 
above,  this  was  not  always  a  clear  choice.  Hence  the  final  choice  was  based  on 
tho  judgement  of  the  computer  operator  with  the  empirical  criterion  of  half  the 
instability  valuo  kept  in  mind.  The  values  cliosen  are  marked  on  Figa.l  to  5, 

2, 2  Frequency  responses 

Tho  fronuoncy  responses  measured  are  defined  ac  the  ratio  of  an  output 
parameter  of  the  system,  say  displacement,  to  an  input  parameter,  say  wind 
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variation,  at  a  particular  frequency*  Since  a  frequency  response  is  only 
valid  for  a  linear  system,  the  convergenco  characteristics  of  the  glide  path 
had  to  be  ignored  and  a  linear  sensitivity  assumed,  For  the  reasons  stated  in 
Section  2.1,  the  sensitivity  chosen  was  that  corresponding  to  a  height  of 
100  ft. 

Three  different  frequency  responses  were  made  for  each  modified  control 
law  and  the  results  obtained  were  compared  with  the  results  from  the  other 
analysis  techniques*  l Measurements  were  taken  to  establish  the  relationship 
between  the  displacement  of  the  aircraft  from  the  glide-path  to:- 

(a)  noise  on  the  I,L,S*  radio  signal, 

(b)  horizontal  wind, 

(c)  vertical  wind 

at  frequencies  in  the  range  0,01  to  1  c/s. 

2.3  Statistical  results 

In  order  to  verify  that  frequency  responses  measured  at  &  fixed  height  were 
representative  of  performance  during  the  approach,  a  set  of  approaches  was 
simulated  for  each  modified  system  in  the  presence  of  random  horizontal  wind 
gusts,  and  the  aircraft  displacement  and  velocity  at  right  angles  to  the  beam 
centre-line  at  the  start  of  ATTITUDE  were  measured  and  analysed*  Fran  these 
measurements,  means  and  standard  deviations  were  calculated  for  each  set  of 
approaches  and  the  results  compared  with  those  of  the  frequency  response 
analysis* 

The  wirri  input  used  in  the  statistical  measurements  was  produced  from  & 
random  whi  te-noise  generator.  The  output  of  the  generator  had  &  Gaussian 
distribution  of  amplitude  with  a  power  spectrum  which  was  sensibly  flat  from 
0.04  cps  to  10  cps.  Since  all  computing  was  done  in  accelerated  time  (time 
constants  decreased  by  a  factor  of  ton)  the  effective  spectrum  was  flat  from 
0.004  cps  to  1  cps* 

The  spectrum  was  shaped  by  passing  the  output  through  a  filter  of  the 

form  ^  — jjy  where  L  is  a  characteristic  turbulence  dimension, 

6 

here  =  1000  ft,  V  is  the  aircraft  airspeed  =  186  ft/sec.  This  gave  a 
e 

filter  of  7 "  ■■■‘v  7*  y:  giving  a  cut-off  at  6  dB  per  octavo  from  0*03  cps.  The 
transfer  function  characteristic  of  the  filter  is  plotted  in  Fig* 9* 
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The  timo  histories  were  intended  to  give  an  overall  qualitative  picture 
of  the  approach  from  beam  capture  to  start  of  ATTITUDE  in  the  presence  of 
wird  shear,  disoretc  gusts  or  in  no  wind*  The  purpose  of  these  time  histories 
was  not  to  differentiate  between  the  systems  but  to  ensure  that  the  other  methods 
of  system  amlysis  had  not  overlooked  important  aspects  of  performaroe. 


MODIFICATIONS  TO  Thd  CONTROL  LAWS 


The  first  retirement  of  any  modification  was  to  improve  the  beam  holding 
in  the  presence  of  tuibulence  without  decreasing  the  stability  at  low  height. 

The  socorri  requirement  was  that  the  modification  should  be  a  simple  one  to 
onginoor.  Sinco  displacement  from  the  beam  is  precoded  by  velocity  and  acceler¬ 
ation,  then  information  about  either  velocity  or  acceleration  normal  to  the 
beam  will  provide  advance  information  on  beam  displacement  errors.  Because  of 
the  consideration  of  engineering  simplicity  it  was  decided  to  investigate  the 
use  of  a  signal  derivod  from  an  accelerometer  positioned  in  the  aircraft  on  the 
centre  of  gravity  so  as  to  provide  a  measure  of  acceleration  along  the  z-axis  of  the 
aircraft.  Velocity  along  the  z-axis  was  obtained  by  integration  of  the  acceler¬ 
ometer  signal.  The  information  derived  from  the  accelerometer  was  used  in 
conjunction  with  that  from  a  pitch  rate  gyro,  to  produce  the  effect  of  displac¬ 
ing  the  accelerometer  from  the  centre  of  gravity.  The  complete  control  equations 
U30d  arc  listed  in  the  Appendix,  but  in  brief  tho  modifications  investigated  were 


the  addition  of: 


(a)  f  (PH) 

(b)  f  (PH  +  P0) 

(o)  f  (DH  + 

(a)  f  (d2h) 

(o)  f  (d2h  +  D2e) 

/  k6  Pv 

to  the  basic  equation  0^  =  -  f  0  +  — ^ — J  • 

Included  in  the  simulation  of  the  accelerometer  was  a  filter  of  time  constant 
0,2  soc,  to  represent  tho  filter  which  would  be  needed  in  practice  to  attenuate 
mechanical  vibration  of  the  accelerometer. 

Tho  possibility  of  adding  a  simple  beam-rate  term  to  the  basic  equation 
was  considered,  but  this  was  rejected  for  a  number  of  reasons.  Fig. 6,  the  plot 
of  the  neutral  stability  bourdary  for  beam-rate  damping  shows  that  a  significant 
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increaso  of  overall  gearing  is  only  possible  if  the  lag  on  the  signal  is  loss 
than  about  0.1  sec.  Unfortunately,  with  such  a  small  lag  the  amount  of  noise 
required  to  cause  the  demanded  pitch  rate  to  exceed  the  limit  of  3  /seo  is  less 
than  12k  micrormrps  at  frequencies  between  1.0  and  2.0  c/s,  see  Pig. 8.  This  is 
less  than  that  .allowed  by  I.C.A.O.  so  that  glide  path  beams  within  the 
specification  would  be  likely  to  cause  rate  limiting.  Further  to  this,  the 
shapo  of  the  stability  boundary  is  so  sharply  peaked  that  even  assuming  no 
limit  problems,  the  practical  difficulty  of  producing  the  required  gearing 
values  would  bo  considerable.  Another  possible  addition  to  the  basic  elation 
that  was  considered  was  that  of  a  simple  pitch-rate  torm,  D0.  However  this 
gave  no  improvement  in  stability  as  is  shown  in  Pig. 8  and  it  was  rejected. 

4  nETMTT.-p  DISCUSSION  OP  TIE  RESULTS 

4.1  Preononcy  responses 

4.1,1  Rosponse  of  the  aircraft  to  boom  signal 

Tho  criterion  against  which  to  compare  the  beam  response  caused  by  the 
various  modifications  to  the  control  law  was  taken  to  be  the  beam  response  of 
the  unmodified  system  using  the  basic  control  law.  The  requirements  for  beam 
rosponse  are  (a)  unity  gain  from  DC  up  to  about  1  o/s  and  (b)  a  sharp  cut-off 
at  about  1  o/s  to  attenuate  the  effects  of  high  frequency  noise.  Since  the 
cloaed-loop  characteristic  equation  of  the  aircraft/autopilot  syst«n  contains 
complex  poles,  such  requirements  oro  not  met  sinco  rosonanoes  produce 
greater  than  unity  gain*’’ Hence  the  requirements  for  the  modified  ays  ton 
are  that  they  increase  bandwidth  without  increasing  the  resonant  gain. 

Pig.  11  to  14  show  the  effects  of  the  various  modifications. 

It  is  clear  that  the  addition  of  D2H  only  or  (D2H  +  D20 )  each  produce 
unacceptable  resonant  gains  at  0.11  c/s  (Pig.  13),  indicating  a  poor  boom  holding 
in  tho  presence  of  disturbances  at  this  froquoncy.  The  addition  of  terns 
DH;  (DH  +  D6)  ard(D^H  +  DH)  all  give  a  wider  bandwidth  with  a  lower  resonance, 
although  tho  presence  of  the  pitch-rate  tom  appears  to  sharpen  the  peak 
slightly. 

4*  1.2  Response  to  horizontal  wind 

Tho  frequency  response  measured  was  tho  amplitude  of  displacement 
produced  by  a  wirci  whose  velocity  varied  sinusoidally,  and  was  plotted  as  feot 
displacement  por  foot  per  seconl  of  horizontal  wind  (Pigs. 15  to  18).  Idoally  it 
is  required  that  the  aircraft  rosponse  to  wind  is  zero  at  all  frocjioncies,  but 
Pig,  15  shows  that  in  the  case  of  the  basic  unmodified  syston  this  is  certainly 
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not  so.  At  gust  frequencies  around  0,1  c/s,  on  amplitude  of  ±5  kts  will 
produoc  deviations  from  glide -path  of  about  ±12  ft  at  a  height  of  100  ft. 

The  spectrum  is  suoh  that  between  periods  of  50  seconds  and  8  seconds,  a 
5  kt  variation  on  the  v/irvi  produces  greator  than  7  ft  deviation  from  the  glide 
slope,  and  since  100  ft  height  is  tho  start  of  ATTITUDE,  this  deviation  will 
carry  through  to  produce  a  touchdown  range  error  of  about  150  ft, 

Tho  modifications  involving  and  (D2H  +  D26),  sec  Pig,  17,  cruise  an 
unacceptable  resonanco  at  0,11  o/s,  although  tho  response  away  from  the  peak 
is  an  improvement  over  that  obtained  with  the  unmodified  control  law  (Fig,  15)# 

A  similar  resonance  was  obtained  on  the  responses  duo  to  beam  signals.  In  the 
case  of  the  modifications  with  DH,  (DH  +  DO)  ana(D^H  +  DH)  (Pigs,  16  and  18)  tho 
curves  have  moved  bodily  to  higher  frequencies  relative  to  the  basic  response 
curve,  and  the  gains  decreased.  This  has  the  effect  of  greatly  improving  the 
responses  at  lower  frequencies  while  giving  a  slight  degradation  at  frequencies 
above  about  0,2  c/s,  i,  o,  about  5  sec  period  gusts.  Pig,  9  shews  that  at 
0,2  c/s  the  amplitude  of  tho  horizontal  wind  is  attenuated  to  about  0,14  that 
at  zero  frequency  so  that  a  slight  degradation  of  performance  at  this  frequency 
is  unimportant,  Tho  greatest  reduction  in  horizontal  wind  response  was  produced 
by  tho  (D2H  +  DH)  modification  (Pig,l8),  tho  reduction  varying  from  a  factor  of 
ten  at  0,003  c/s  to  a  factor  of  5  at  0,1  c/s.  Above  this  frequency  the  improve¬ 
ment  becomes  less  until  at  0,2  c/s  the  response  curve  of  tho  basic  system  crosses 
that  of  tho  modified  control  lav/, 

4,1, 3  Response  to  vortical  wind 

The  response  of  the  system  to  vertical  winds  (^igs.^  to  22)  was  measured 

as  a  function  of  tho  height  error  caused  by  variation  of  tho  vortical  wind 

speod.  The  results  as  given  by  the  vertical  wind  responses  are  identical  to 

2  2  2 

those  of  the  horizontal  wind  case  in  that  D  H  and  (j)  H  +  D  6)  give 

2 

exaggerated  maxima  while  the  least  response  is  from  (D  H  +  DH),  The  latter 

modification  gives  a  reduction  factor  of  about  2,7  at  frequencies  less  than 

about  0,15  c/s  when  its  response  curve  again  converges  with  and  crosses  that  of 

the  basic  system,  7/ith  the  basic  system  tho  maximum  displacement  from  glide 

2 

path  is  2,8  ft/ft/sec  for  0,1  o/s  gusts,  Tho  greatest  displacement  v/ith  (D  H+DH) 

occurs  at  a  frequency  of  0,1 85  c/s  when  the  displacement  is  1,3  ft/ft/sec.  The 
2  2  2 

modifications  with  D  H  and  (D  H  +  D  0)  give  a  degraded  response,  the  displacement 

2  2 

produced  by  1  ft/sec  gu3ts  being  4  ft  for  (D  H  +  D  0)  at  0,11  c/s. 
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1^2  Statistical  treatment  of  response  to  random  horizontal  wind. 

In  order  to  verify  the  predictions  given  by  the  frequency  re^jonaee, 
simulated  approaches  wore  made  for  each  modification  state  in  the  presence  of 
raniom  horizontal  wind  with  an  rms  velocity  of  4  ft/soc  equivalent  tp  a  moan 
winl  speed  of  about  12  kts  and  a  simple  statistical  analysis  made  of  the 
displacement  and  velocities  perpendicular  to  the  glide  slope  at  tho  start  of 
ATTITUDE.  Twonty-five  simulated  approaches  were  made  with  each  control  law, 
and  each  approach  was  stopped  at  a  height  of  100  ft  and  the  above  parameters 
recorded  from  a  digital  voltmeter.  Ihe  means  and  standard  deviations  of  tho 
displacements  and  velocities  were  than  calculated,  and  from  these  an  ostimato 
of  touchdown  range  error  duo  to  random  wind  variations  was  made* 

Previous  work2  has  shown  that  the  touchdown  rang#  error  produced  by  a 
vertical  velocity  error  of  1  ft/sec  at  the  start  of  ATTITUDE  is  about  175  ft 
for  a  nominal  rate  of  descent  of  10  ft/sec  (a  3°  glide  path  and  an  approach 
speed  of  186  ft/sec  gives  a  nominal  rate  of  9.7^  ft/sec).  Assuming  such  a 
range  error  due  to  velocity  error,  and  assuming  that  displacement  errors  ore 
carried  through  geometrically  to  give  19»1  ft/ft  of  beam  error  for  a  3  glide 
slope  Table  1  gives  the  Stanford  deviation  of  touchdown  range  caused  by 
horizontal  wild  for  the  different  control ‘laws,  and  Fig.  23  gives  examples  of  the 
approaches  in  the  presence  of  random  horizontal  wind  of  rms  velocity  4  ft/seo* 

Table  1 


Additional 

control 

function 

Velocity 
error  (Dh) 

(ft/sec) 

Displacement 
error  (h) 

(n) 

Touchdown  error 
in  ft  duo  to 

Total 

touchdown 

error 

(«) 

Dh 

h 

None 

a.  94 

3.44 

1 6A 

06 

176 

DH 

0.54 

94 

29 

98 

(DH+rfi) 

0.39 

67 

18 

70 

(DH+  dV; 

0.41 

72 

21 

mm 

d2h 

0.75 

132 

84 

HI 

(D2H+D20) 

0.71 

124 

i 

45 

1 

All  the  above  numbers  represent  standard  deviations. 


Fig,  2d  illustrates  in  diagrammatic  form  tho  difference  in  range  error 
produced  by  the  different  modifications. 


12 


i*3  Time  history 

Simulatod  beam  captures  and  approaches  were  made  from  a  height  of 
2000  ft  for  each  control  law  investigated  in  threo  different  sets  of  wind 
conditions 

(i)  sero  wind, 

(ii)  tailwind  shear  (i.e.  tailwind  whose  speed  is  dropping  with  height), 

(ill)  headwird  shear  with  a  horizontal  gust  of  5  ft/sec  at  300  ft  height* 

In  each  instance  it  was  considered  that  if  the  beam  joining  manoeuvre 
for  a  particular  control  law  was  not  grossly  different  from  that  of  the  others, 
then  small  differences  in  behaviour  were  only  important  at  close  range,  i.e. 
at  heights  below  500  ft, 

1*3.1  Approaches  in  zero  wind 

As  seen  from  fig. 25  all  of  the  modifications  to  the  control  law  produce 
improvements  in  the  beam  joining  in  zero  wind  conditions.  Normal  aircraft 
instruments  would  indicate  that  the  aircraft  was  established  if  the  error  from 
the  beam  was  less  than  about  15  microamps;  this  would  mean  that  all  of  the 
systems  were  cstablisned  iimediately.  Hence  all  of  the  modifications  produced 
satisfactory  beam  joining  in  sero  wind  conditions. 

4. 3, 2  Ar  ^roaches  in  the  presence  of  a  tailwind  shear 

The  wind  simulation  was  a  tailwind  of  50  ft/scc  at  a  height  of  2000  ft 
falling  linearily  to  20  ft/soc  at  ground  level,  and  Pig. 26  shows  that  if  the 
criterion  of  15  microamps  is  taken,  then  those  systems  using  cither  no  extra 
tsm  of  DH  dancing  only  are  not  established  until  below  500  ft  height.  The 
other  ays  toms  all  roach  the  15  microamp  criterion  by  a  height  of  about  1200  ft 
although  the  system  using(DH  +  D6)  damping  only  just  remains  inside  the 
15  micronmp  limit  after  the  second  crossing  of  the  boam  centre-line. 

4,3#3  Approaches  in  the  presence  of  a  headwind  shear  and  a  gust 

Pig,27  shows  the  combined  affect  of  a  wind  shear  and  a  horizontal  wind 
gust  from  ahead,  the  gust  being  a  step  of  5  ft/sec  injected  at  a  height  of 
300  ft.  The  wind  Simula  tod  for  the  headwind  shear  was  a  headwind  of  50  ft/sec 
at  a  height  of  2000  ft,  falling  linearly  to  20  ft/sec  at'  ground  level. 

All  of  the  systems  achieved  good  beam  joining  performances,  even  the 
unmodified  system  resulting  ir  the  glide  path  error  signal  being  less  than 
15  microoirps  by  a  height  of  1400  ft.  Those  control  systems  using  (D  H+D  6), 
(D2H  +  DH)  and  D2H  only  were  within  this  limit  Immediately  and  kept  to  within 
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5  miorownps  after  the  third  beam  crossing.  This  is  much  better  performance 
than  that  achieved  in  a  tailwind  shear,  but  it  should  be  noted  that  the  change 
in  performance  is  entirely  due  to  the  sign  of  the  wind  direction  at  2000  ft 
height  rather  than  to  the  sign  of  the  shear  itself.  Even  in  zero  wind  the 
aircraft  overshoots  the  glide  path  and  a  tailwind  makes  the  overshoot  worse, 
whilst  a  headwind  reduces  it.  Nevertheless  since  the  general  tondenoy  is  for 
wind  strength  to  decrease  as  height  is  reduced,  it  is  considered  that  the 
method  of  simulating  wind  shear  represents  a  disturbance  likely  to  be  encountered 
in  practice. 

The  response  of  the  systems  to  a  step  gust  followed  approximately  the 

results  prodicted  by  the  frequency  responses  except  for  the  response  of  the 
2  2  2 

systems  using  (D  H  +  D  0)  and  0  H  damping,  vrtiose  gust  responses  were 
predicted  to  be  worse  than  those  of  the  basic  system.  It  must  be  realised, 
however,  that  the  gust  was  injected  at  a  height  of  300  ft  to  allcw  time  for  a 
complete  recovery  while  the  frequency  responses  were  measured  at  a  height  of 
100  ft.  Frequency  responses  measured  at  heights  greater  than  100  ft  show  that 
the  vory  sharp  resonance  in  the  responses  of  the  system  using  D^H  and  (D2H  +  D2e) 
falls  off  rapidly  with  increase  in  height  (see  Fig,l6)  and  frequency  responses 
measured  at  300  ft  agree  vdth  the  results  of  the  gust  injection  at  300  ft, 

4,4  Effect  of  an  error  on  the  datum  rate-of-descont  in  the  control  laws 

using  DH 

The  use  of  a  rate-of -descent  term  for  damping  in  the  control  equation 
requires  that  the  velocity  signal  derived  from  the  accelerometer  bo  biased, 
so  that  zero  eignal  is  obtained  when  the  aircraft  is  descending  at  a  rate 
appropriate  to  the  glide -path  angle  and  approach  speed  of  the  aircraft. 

Since  beam  angles  will  vary  between  about  2,5°  and  3*5°  and  aircraft 
ground  spood  will  vary  according  to  wind  conditions,  approaches  wore  recorded 
for  the  control  law  modifications  involving  DH,  (DH  +  DH)  aid  (DH  +  DO)  in 
various  wind  conditions,  with  errors  on  the  datum  rate  of  descent  of  about 

±30' a,  equivalent  to  beam  angles  of  3°  ±  1°.  This  is  also  oquivalent  to  an 

actual  beam  angle  of  3°  ±  0.5°  plus  a  head  or  tailwird  of  28  ft/seo  (l6,5  kts). 

The  sign  convention  used  is  as  follows:-  if  the  datum  rate  of  descent 
set  into  the  pitch  computer  is  less  than  that  roquirod  for  the  glide-path 
angle,  then  the  error  is  defined  as  negative.  Hence  an  error  of  -3 O/o  on  a 
glide  path  of  nominal  rate  of  descent  9,73  ft/sec  means  that  the  datum  is  set 

for  a  rate  of  do  scent  of  6.82  ft/sec.  Figs.  28  to  30  show  the  effect,  for  a 

glide-path  with  a  nominal  rato  of  descent  of  9,75  ft/sec,  of  an  error  of 
t  3  ft/sec. 
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Fig. 28  shows  the  effects  of  such  errors  in  the  absence  of  any  outside 
disturbance,  other  than  the  beam  join  manoeuvre,  and  comparison  with  Figs.25 
arri  26  indicate  the  similarity  of  the  effect  of  wind  shears  from  head  and  tail 
with  no  datum  error.  This  is  only  to  be  expected  since  both  effectively 
produce  ground  speed  errors  if  the  airspeed  is  assumd  to  be  constant.  When 
a  30?$  error  on  the  datum  rate-of -descent  is  coupled  with  either  a  head  or 
tailwird  shear,  Figs. 29  and  30,  the  effect  of  the  error  is  added  algebraically 
to  the  shear  so  that  either  an  improvement  or  a  degradation  in  the  beam  join 
manoeuvre  may  be  achieved.  In  both  cases  the  effect  of  the  error  is  only 
evident  during  the  first  part  of  the  approach,  there  being  little  difference 
in  the  path  flown  with  perfect  datum  ard  with  a  datum  error  below  a  height 
of  about  600  ft. 

Even  with  as  large  an  error  as  30#  on  the  datum,  the  performance  in 
glide-path  capture  using  DH,  (DH  +  DO)  or  (D2H  +  DH)  in  the  control  law  is 
comparable  to  that  of  the  basic  control  law  as  used  at  present  in  the  Varsity/ 
system. 

5  CONCLUSIONS 

The  purpose  of  this  study  was  to  investigate  the  possibility  of  improving 
the  aircraft’s  performance  on  the  I.L.S.  glide  path  and  so  reduce  the  effect 
of  wind  variations  on  the  scatter  of  the  touchdown  point.  The  results  obtained 
confirm  that  such  an  improvement  is  possible  for  the  simulated  system  using 
the  additional  terms  (DH  +  D0)  or  (D2H  +  DH)  in  the  control  law. 

Yfhile  those  terms  do  not  produce  the  best  beam  capture  in  the  presence 
of  winds  hears,  it  is  considered  that  the  greater  improvement  in  the  presence 
of  v/ind  gusts,  coupled  to  the  higher  probability  of  encountering  wind  gusts 
than  that  of  encountering  wind  shears,  recomnends  the  use  of  either  (DH  +  W) 
or  (D2H  +  DH). 

The  use  of  either  of  these  terms  has  been  shown  to  reduce  the  standard 
deviation  of  touchdown  range  due  to  horizontal  vand  turbulence  in  the  ratio 
3:2  and  even  allowing  for  practical  deficiencies,  a  considerable  reduction 
should  be  achievable. 

A  flight  programme  is  now  in  hand  to  investigate  the  practical  validity 


of  the  conclusions. 
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Appendix 

GKNERAL  EQUATIONS  USED  IN  THE  S  EMULATION 
Control  law  modifications 

The  basic  control  law  as  used  at  present  in  the  B.L.E.U.  system  iss- 

-K, 


6c  "  (1  7  0.2  D)Ji  7  0.5  d)  [P  +  K6  d] 

The  modifications  to  the  control  law  produced:- 


-s  ,  1 

rp  +  K  i  +  "ioi  “  1 

'(1  +  0.2  DM1  +  0.5  d;  I 

+  K6D  +  (1  +  0.2  dIJ 

-K<  ,  I 

Pa  ,  K  2.  +  K191  .  .a  [dh  ♦  K  D0]*] 

(1  +  0.2  D)(1  +  0.5  D)  1 

l_P  +  o  D  +  (1  +  0.2  d)  1  105  JJ 

_K?  —J 

fp  +  K  £+  K102  D2j_] 

(1  +  0.2  D)(l  +  0.5  D)  1 

+  1  D  (1  ♦  0.2  D) J 

~S 

(1  7  0.2  D)(l  +  0.5  d) 

[p  +  K6  D  +  (1  +  0.2  D)  LV  n  +  103  JJ 

-K5 

fp  +  K  £  ♦  K101  PH  *  K102  D2H1  . 

(1  ~  0.2  DKi  +  0.5  d) 

[P  +  K6  D  +  (l  7  0. 2  D]  J 

Kg  was  fixed  at  l/30  seo”1  for  all  of  the  equations. 

f  the  other  coefficients 

,  optimised  for  each  equation,  were:- 

K5  =  0.02%* 

a  0,03  /M-®-,  “^01 

=  7  pa/ft/sec 

K5  =  0.04%a,Kl01 

=  7  Ha/ft/sec,  =  1.0  ft/seo/%ec 

K5  =  0.03%a,  KiQ2 

2 

*  3  ^a/ft/sec 

K5  =  0.04%a,K102 

=  5  na/ft/seo2,  K1QJ  a  0.175  ft/se<^°/eec? 

K5  =  0.05%a,  KiQ1 

=  7  na/ft/aec  ,  KlQ2  =  3  na/ft/seo2 
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Appendix 


Aerodynjunlg  equations 

After  aub.titution  for  the  aerodynamic  derivativea  of  the  Varsity 
in  the  approach  configuration,  the  equations  are:- 

«*.««  (7)  •  0.O2K  («•«.)-  °-”8  *  *  '  - 

d«  de+ol234T,  =  0 

Equation  (8)  0.1068  (u  ♦  \)  *  0.938  «  +  dt  ♦  dt  dt 

.  .  n  ni  *  2  2  «  +  ♦  1.481  dT  +  524  *1  =  0 

Equation  (9)  0.474  dt  +  2,2  7?  dt 

Enuatioreof  motion 

V  (0-“) 

.  .  dH  IsJ - 1 

Equation  (10)  dt  =  57.3 

(V  -  W)e 

dh  _  dH  .  l_e - 

Equation  (11)  Jt  ”  dt  57.3 

CauMml  f  123  B  =  R  -  [V  -  *  ♦  (“  +  %)]  t 


Equation  (12)  R  -  R0  -  ^  “ 

a  Rp 

Equation  (l3)  0  3  R  h 

VH  nd  equations 

H 

Equation  (14)  Shear  :  VI  =  »2  ♦  (»,  -  W2)  H" 

Equation  (15)  Longitudinal  wind  :  u^  =  W  '  W1 

We 

Equation  (If)  Vertical  wind  :  «„  =  v# 


n^r.tTf.1  aquations 

G  (i  +  °-3  -  6J  *  G2  1 

Equation  (17)  Basic  autopilot  .  u  '  G1  (l  +  Q>1  D)2  L  8 

Equation  (l8)  P  =  £(9  “  V  +  ^dt] 

-K^  (p  +  Kg  g  +  f(DH,D2H,De,D28)j 
Equation  (19)  Locking  unit  :  P0  -  (l  +  0.2  D)(1  7  0.5  u) 


Equation 


Equation  (20)  Autothrottle  :  -  = 


-T  r  (u  +  uwh 

7  OTw  KTdt[1u  ‘  *  t7  d  J 


Appendix 
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Values  of  parameters 

e  =  +3° 

R  =  38200  ft 

o 

W  =  50  ft/sec 

'.Yg  =  20  ft/sec 

V  =  186  ft/seo 

o 

II  =  2000  ft 

o 

sRp  =  18000  millivolts 

=  2.0  degrees /degree 

G2  a  1/15  sec"1 

=0.1  (equivalent  to  about  180  lb  At  for  Varsity  of  mass  1055  slugs) 
-r  n.05  seo"’ 

Tj  =  0.35  (=  590  lb/degree) 
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SYMBOLS 

0  *  commnded  pitoh  angle  in  degrees 

o 

u  m  forward  airspeed  change  in  ft/sec 


a  n  in  degrees,  where:- 

w  s  aircraft  airspeed  along  its  z  axis 
Ve  s  da  tun  airspeed 

6  s  pitoh  angle  in  degrees,  relative  to  the  earth's  horizontal 

Tj  s  elevator  angle  in  degrees 

H  =  height  of  aircraft  above  ground 

E  =  initial  height  of  aircraft  above  ground 
o 

h  .  perpendicular  height  of  aircraft  above  or  belwr  the  glide  slope 

R  s  range  from  glide  path  aerial 

R  s  starting  range  of  aircraft  from  glide  path  aerial 
4  o 

W  =  wind  speed  at  a  height  H 
3  wind  speed  at  a  height  Hq 

\72  a  wind  speed  at  ground  level 

cVG2*K5,K6  ***  *utoPilot  Parametera 


Wt3 


are  auto  throttle  parameters 
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FIG.  5  NEUTRAL  STABILITY  BOUNDARY  FOR 
(D2H+D2e)  DAMPING 
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FIG.  6  NEUTRAL  STABILITY  BOUNDARY  FOR 
D0  DAMPING 
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Fig. 7 


FIG.  7  NEUTRAL  STABILITY  BOUNDARY  FOR 
D0  DAMPING 
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FIG.8  NOISE  INPUT  REQUIRED  TO  EXCEED  THE  3°/SEC  PITCH  RATE  LIMIT  WITH  D0  DAMPING 
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FIG. 9  CHARACTERISTIC  OF  THE  SIMPLE  FILTER 
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FIG. II  FREQUENCY  RESPONSE  TO  RADIO  SIGNAL 
OF  THE  BASIC  CONTROL  SYSTEM 
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FIG.  12  FREQUENCY  RESPONSE  TO  RADIO  SIGNAL 
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FIG.I3  FREQUENCY  RESPONSE  TO  RADIO  SIGNAL 
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FIG.  14  FREQUENCY  RESPONSE  TO  RADIO  SIGNAL 
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FIG.  16  FREQUENCY  RESPONSE  TO  HORIZONTAL 
WIND  VARIATIONS 
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FIG.  17  FREQUENCY  RESPONSE  TO  HORIZONTAL 
WIND  VARIATIONS 
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Fig. 19 


FIG.  19  FREQUENCY  RESPONSE  TO  VERTICAL  WIND  OF 
THE  BASIC  CONTROL  SYSTEM 
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FIG.  20  FREQUENCY  RESPONSE  TO  VERTICAL 
WIND  VARIATIONS 
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FIG.  21  FREQUENCY  RESPONSE  TO  VERTICAL 


WIND  VARIATIONS 


Fig. 22 
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FIG  22  FREQUENCY  RESPONSE  TO  VERTICAL 
WIND  VARIATIONS 


Fig. 24 
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FIG  24  EFFECT  OF  CONTROL  LAW  MODIFICATIONS  ON 
STANDARD  DEVIATION  OF  TOUCHDOWN  RANGE 
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F!6.25  EFFECT  OF  ADDITIONAL  DAMPING  TERMS  ON  THE  GLIOEPATH 


PERFORMANCE  IN  ZERO  WIND. 
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FIG.  26  EFFECT  OF  ADDITIONAL  DAMPING  TERMS  ON  THE  GLIDEPATH 
PERFORMANCE  IN  THE  PRESENCE  OF  TAIL  WIND  SHEAR 
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FIG.27  EFFECT  OF  ADDITIONAL  DAMPING  TERMS  ON  THE  GLIDEPATH 
PERFORMANCE  IN  THE  PRESENCE  OF  HEAD  WIND  SHEAR 
AND  A  GUST  OF  5  FT/SEC  AT  300  FT 
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FIG  28  EFFECT  OF  ERROR  ON  RATE  OF  DESCENT  OATUM  FOR  A 
3°  GLIDEPLATH  WITH  NO  EXTERNAL  DISTURBANCE. 
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Fig. 29 
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FI6.29  EFFECT  OF  ERROR  ON  RATE  OF  DESCENT  DATUM  FOR 
A  3*  GUDEPATH  IN  THE  PRESENCE  OF  TAILWIND  SHEAR 
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Fig. 30 


FOR  A  3°  GLIDE  PATH  IN  THE  PRESENCE  OF  HEADWIND  SHEAR 


